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We construct new tests of perturbative QCD by considering a hypothetical r lepton of arbitrary 
mass, which decays hadronically through the electromagnetic current. We can explicitly compute 
its hadronic width ratio directly as an integral over the e + e~ annihilation cross section ratio, R e + e - ■ 
Furthermore, we can design a set of commensurate scale relations and perturbative QCD tests by 
varying the weight function away from the form associated with the V — A decay of the physical r. 
This method allows the wide range of the R e + e - data to be used as a probe of perturbative QCD. 



The hadronic width of the r lepton is potentially one 
of the most important sources for the high precision de- 
termination of the coupling ajjg of QCD [|],|| . The per- 
turbative QCD (PQCD) analysis of the r width has been 
refined by constructing moments of hadronic decay distri- 
butions which minimize sensitivity to the low energy part 
of the hadronic spectrum g]. However, it is still uncertain 
whether the r mass is sufficiently high to trust PQCD, 
particularly due to the strong distortion of hadronic final 
state interactions . 

In this paper we construct new renormalization 
scheme-independent tests of PQCD which we can apply, 
not only to the physical r lepton, but also to a hypo- 
thetical r lepton of arbitrary mass which decays hadron- 
ically through the vector current. Such hypothetical r 
leptons, with masses M < M T , have already been con- 
sidered in ref. || . We can obtain empirical values for the 
hypothetical lepton's hadronic width and moments di- 
rectly as integrals over the measured R e + e - = cr(e + e~ — > 



hadrons)/cr(e H 



). As we shall show, these 



tests are fundamental properties of QCD which can serve 
as necessary conditions for the applicability of perturba- 
tion theory. 

Quantum field theoretic predictions which relate phys- 
ical observables cannot depend on theoretical conven- 
tions such as the choice of renormalization scheme or 
scale. The most well-known example is the "generalized 
Crewther relation" || in which the leading twist PQCD 
corrections to the Bjorken sum rule at a given lepton 
momentum transfer Q 2 are inverse to the QCD correc- 
tions to R e + e - at a corresponding CM energy squared, 
s* = s*(Q 2 ), independent of renormalization scheme. 
The ratio of the scales s* /Q 2 has been computed to NLO 
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in PQCD. Such leading-twist predictions between observ- 
ables are called "commensurate scale relations" and are 
identical for conformal and nonconformal theories Q . 

Another important example is the commensurate scale 
relation between the PQCD correction to the r lep- 
ton's width ratio, R T — L(t~ — > v T + hadrons)/r(r _ — > 
v T e~v e ), and those to R e + e - ■ Assuming for now / mass- 
less flavors, PQCD yields 



R P 



f 
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where a R can be written as a series in a s /ir in some renor- 
malization scheme. Note that a R is an effective charge 
because it satisfies the Gell-Mann-Low renormaliza- 
tion group equation with the same coefficients /3o and (3\ 
as the usual coupling a s (differing only through the third 
and higher coefficients of the /3- function) . Similarly we 
can define an effective charge a T as follows 



R T (M T ) = R° T (M T ) 



1 



a T (M T ) 



Leading-twist QCD predicts 

a T {M T ) = a R (Vs*) 



(2) 



(3) 



to all orders in perturbation theory. The ratio of the 
commensurate scales is known in NLO PQCD: 
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This result was originally obtained in |6| by using NNLO 
predictions for a R and a T obtained in the MS scheme 
and eliminating ctjjg. However, as we shall show here, 
the QCD prediction for \fs* jM r also follows from the 
fact that both effective charges evolve with universal /3o 
and (3\ coefficients. The fact that R T can be expressed 
asQ 
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1, as in 



Note that in 



order to include NNLO corrections in a T , we must modify the 
O(o4) coefficient of as by setting (J~) f qf) 2 = 0. 
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implies, by the mean value theorem, that a R and a r are 
related by a scale shift. However, the prediction for the 
ratio y/s*/M T in eq.([IJ) is a specific property of PQCD. 

A definitive empirical test of the commensurate rela- 
tion, eq.(0), is problematic since there is only one t lepton 
in nature, and its mass seems uncomfortably low for tests 
of leading-twist QCD. However, we can construct new 
tests of PQCD by considering a hypothetical r lepton 
of arbitrary mass M which decays hadronically through 
the vector current. Then we can explicitly compute its 
hadronic width ratio as an integral over the measured 
R e + e -. Furthermore, we can design a set of commen- 
surate scale relations and PQCD tests by varying the 
weight function away from the form associated with the 
V — A decays of the physical r. Thus we can use the 
full range of the R„ novel test of PQCD. As we 

shall show, such a test must also take into account spe- 
cific effects attributable to the ss, cc, bb quark thresholds. 
Also, following |l we shall smear the annihilation data 
in energy in order to eliminate resonances and other dis- 
tortions of final state interactions. By smearing R e + e - 
over a range of energy, AE, we focus the physics to the 
time At = 1/AE where an analysis in terms of PQCD 
quark and gluon subprocesses is appropriate. Therefore, 
this method can also be interpreted as an additional test 
of duality. Scheme-independent relations between R e + e - 
and r decay have also been recently discussed in M . 

Given a R , we can construct new effective charges as 
follows: 
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We can choose f(x) to be any smooth, integrable func- 
tion of x = s/M 2 . ( For the particular choice, f T (x) = 
(1- x)' (l + 2x), a R is simply a T .) The mean value the- 
orem then implies 



a f {M) = a R (Js* f ), 



< s* f < M . (7) 



Dimensional analysis ensures that y/~Sj = A/ ill, where 
A/ possibly depends on a R . To obtain an estimate for 
Xf we consider the running of a R up to third order: 



<*r(Vs) = a R (M) _ fa . fj_\ ( a R {MY 
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We substitute for a R in eq. (§) to find 



a/(M) = a R (M) _ fa ( h\ fa R (M) 

TT TT 4 W 



+ 



1 

16 



1% 



V 77 
h 

In 



a R (M) 



(9) 



where /; = J Q f (x) (In x) 1 d x . By setting s = s* in eq. (||) 
and comparing with eq. (^|), we extract 



Xf = exp 
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Note that if /(x) is positive on the interval [0, 1], then 
h/Io is negative as expected. Using f T (x), eq. ( |l(i| ) 
is nothing but eq. (Q). Also, since A/ is a constant to 
leading order, af should satisfy the same RG equation as 
a R with the same coefficients (5q and In other words, 
ay is an effective charge. 

We can now study integrals over R e + e - data with dif- 
ferent weight functions f(x) and vary M to see whether 
we obtain the PQCD behavior. In general, the weight 
function f(x) should be chosen to suppress the low en- 
ergy region, where non-perturbative effects are impor- 
tant. Thus, in the following, we will set f(x) = x k , where 
k is some positive number. In such a case, we have that 



(M) = a R (X k M) with A fe = e W2/ 0fc; (n) 



where I\ k = J^x k \nxdx and Iq^ = x k dx. Note 
that as k increases, Iik/^Iok — * 0, and, therefore, \fs* — > 
M . For very large k, we lose sensitivity to the details of 
PQCD. It is particularly interesting to use such a test to 
probe the energy region close to the r mass M T . 

The main difficulty in comparing with R e + e - data is 
that we can no longer consider massless flavors and that 
we observe hadrons instead of quarks. 

Following U the effect of quark masses can be approx- 
imately taken into account if we use: 
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where v% = yl — 4m| /s is the velocity of the initial 
quarks in their CM frame. The t?j(3 — vf)/2 factor is 
the parton model mass dependence and g(v) is a QCD 
modification |10| of the Schwinger positronium correc- 
tions [ pT[ . In principle, all these corrections spoil the 
relation in eq. (|TT|) . However these factors are unity for 
energies well above their corresponding thresholds. 

Nevertheless, we still cannot compare directly with 
the data since there is no direct correspondence between 
quark and hadronic thresholds. To obtain a meaningful 
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comparison we have to smear both the PQCD results and 
the data. Following || we define smeared quantities as 
follows: 



R(V») = - 



(s - s') 2 - A 2 



(14) 



Note that in the A — > limit, we recover the original 
quantity. In what follows we use the standard value 
A = 3 GeV 2 [p|p^] . The smearing effect can be seen 
by comparing Fig.l, which shows an interpolation of the 
R e + e - data, |15|| , with Fig. 2. For completeness, we also 
include in Fig. 2 the smeared results from NLO PQCD 
and from the naive parton model (a.R = 0). 
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Figure 1: Interpolation of the central values of R e + e - data 
|l3| . Narrow resonances are taken into account using their Breit- 
Wigner form. Note that there seems to be a discrepancy in the 
central values of experiments between 5 and 10 GeV, that above 
20 GeV we have two or three clearly different central values at the 
same \fs, and that the point at 13 GeV is much higher than other 
nearby data. 



bottom 



A=3 GeV" 




Re* e - from data 
IV e - from NLO QCD 
Ro (parton model) 



1.2 



10 



20 



35 



Vs (GeV) 

Figure 2: Effect of smearing on _R e + e _ . 

In order to integrate over R e + e - , we need to interpo- 
late, but not fit, the data. Note that any fit using the 
QCD functional dependence will always satisfy the com- 



mensurate scale relations, even if its quality is poor. To 
avoid this bias, we have interpolated the central values of 
the data by means of "r-term simple moving averages" 
up to 30 GeV (to avoid electroweak contributions). That 
is, if we have a series of raw data z\, ...z n , we obtain the 
new set of smoothed data Y^j=o w j Zt ~i ^ or * = r i--- n i 
with Y^j—Q Wj — 1. We have used r ranging from 2 to 6 
for different energy regions and our moving averages are 
"simple" because all the weights Wj are equal. Finally, 
the resulting smoothed data have been interpolated us- 
ing cubic-splines. In addition, the narrow resonances that 
do not appear in Fig.l are implemented using the Breit- 
Wigner formula. 

We have thus eliminated the QCD biases up to 30 GeV. 
Above that energy we have matched a logarithmic func- 
tion whose functional dependence is inspired by QCD, 
but its contribution in the smearing integrals is negligi- 
ble for small a/s. 

Unfortunately, we cannot extract directly the effective 
charges from their corresponding smeared ratios since 
they are multiplied by other functions inside the smear- 
ing integral. However, using eqs.(p^) and (|lj), we define 
smeared charges: 



Re+e-(Vs) - Rpjy/s) 

Rsch(Vs) 



(15) 



and similarly for otk- In the massless A — > limit we 
recover the standard effective charges. We expect the 
smeared charges to satisfy eq. (^) in energy regions where 
the threshold corrections can be neglected. 
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Figure 3: Comparison between an(\/~s*) and different ctfc mo- 
ments at M = \fs* The dotted line shows how the agreement 
is spoilt if we do not shift y/s* to M. 



In Fig. 3 we compare cxr at ys* with ctk moments at 
M = \/s*/Afc. For ao the agreement is poor, since the 
low energy region is not suppressed enough. But for ct\ 
we find a reasonable agreement in several regions, and 
we also show how this agreement disappears if we do 
not shift the argument of ct\ from ys* to M — s/s*/\i. 
Starting from higher energies, we find above 30 GeV that 
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commensurate scale relations are satisfied almost identi- 
cally, which is not surprising since above that energy we 
have fitted with a QCD inspired behavior. From f 5 GeV 
up to 30 GeV different experiments have measured rather 
different central values at very similar, or even the same, 
energies. The smooth interpolation of these points pro- 
duces artificial oscillations around the mean values of the 
data. As far as these oscillations are centered on the at 
curves, there is a reasonable agreement, given the quality 
of the data. In the region between 5 and 10 GeV there 
seems to be some controversy about the compatibility be- 
tween different experiments (see Fig.l and ref. H])- It 
has become standard not to use the older data (which is 
higher) as we have done in Fig.l. Although the more re- 
cent data may be compatible within their experimental 
errors with the QCD expectations, their central values 
are systematically lower, which is why eq. ( |ll|) does not 
seem to hold. Once there are more accurate data, the 
tests we are proposing, together with a thorough error 
analysis, will shed light on this situation. 



imprecision of much of the R e + e - data. Smearing the 
data over an energy range helps but does not totally re- 
move the effects due to final-state interactions. Quark 
threshold distortions are partially alleviated by using the 
Schwinger corrections at small velocity, but the domain 
of non-relativistic velocity introduces its own complica- 
tions, including sensitivity of the running coupling to the 
soft a s m q scale |l5j . Remarkably, the mass range of the 
physical r lepton is potentially clear of the finite quark 
mass effect since it is well below the cc threshold. How- 
ever, it is clear that higher precision measurements of 
R e + e - throughout the energy domain below the Z° bo- 
son are needed. 

We are indebted to M. Swartz for assistance with the 
experimental data. J.R.P. thanks the Spanish Ministerio 
de Education y Cultura for financial support. 
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Figure 4: Comparison between &r(-\/s*) and different <5j. mo- 
ments at M = \fs* /Afc in the low energy region. 

The low energy region is shown in Fig. 4 in more detail. 
Taking into account that we are only using LO QCD and 
central data values, the agreement between the shaded 
regions looks quite satisfactory. This is encouraging for 
the real t lepton, which sits in a region where PQCD 
results may be applicable since it is primarily sensitive 
to the light u,d,s flavors. Nevertheless, by looking at 
energies ^/s ~ 1.5 GeV, our results seem to support the 
claims that the R e + e - data could be 6-7% lower than the 
QCD expectations in that region 

The commensurate scale relations connecting the mo- 
ments of the lepton hadronic decay spectrum to R e + e - 
derived here are basic scheme- independent tests of 
PQCD, depending only on the the universal terms of the 
[3 function. We have seen, however, that a direct compar- 
ison with data is problematic because of several factors 
such as the distortions of narrow and broad resonances, 
the physical effects of the quark pair thresholds and the 
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